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Qh Discrete systems with long-range interactions are considered. Continuous medium 

models as continuous limit of discrete chain system are defined. Long-range interactions 
of chain elements that give the fractional equations for the medium model are discussed. 
The chain equations of motion with long-range interaction are mapped into the contin- 

^ . uum equation with the Riesz fractional derivative. We formulate the consistent definition 

of continuous limit for the systems with long-range interactions. In this paper, we con- 
sider a wide class of long-range interactions that give fractional medium equations in the 
continuous limit. The power-law interaction is a special case of this class. 
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1 Introduction 

Derivatives or integrals of noninteger order [H El El HI [5] have found many applications in 
recent studies in mechanics and physics [6j [7J El El [H)] . Equations which involve derivatives or 
integrals of non-integer order are very successful in describing anomalous kinetics and transport 
and continuous time random walks [HI [121 131 ELU E5] . Usually, the fractional equations for 
dynamics or kinetics appear as some phenomenological models. Recently, a method to obtain 
fractional analogues of equations of motion was considered for sets of coupled particles with a 
long-range interaction [hoi (TTJ dHJ [TJJ]. Examples of systems with interacting oscillators, spins 
or waves are used for many applications in physics, chemistry, biology [2U1 ED E21 E21 EH ESI 
EH E3 EH EH EH1 EH E21 1331 El- In the continuous limit, the equations of motion for discrete 
systems give the continuous medium equation. The procedure has already been used to derive 
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fractional sine-Gordon and fractional wave Hilbert equation [TBI HE]) to study synchronization 
of coupled oscillators [17] , to derive fractional Ginzburg-Landau equation [T7J and for chaos in 
discrete nonlinear Schrodinger equation [15] . In [16], [181 EE3 US] > om y the power-law long-range 
interactions are considered. In this paper, we consider a wide class of long-range interactions 
that give fractional medium equations in continuous limit. The power-law interaction is a 
special case of this class. 

Long-range interaction (LRI) has been the subject of investigations for a long time. An in- 
finite one-dimensional Ising model with LRI was considered by Dyson [20] . The <i-dimensional 
classical Heisenberg model with long-range interaction is described in [211 122] , and their quan- 
tum generalization can be found in [23] [23]. Solitons in a one- dimensional lattice with the long- 
range Lennard- Jones-type interaction were considered in [2B]. Kinks in the Frenkel-Kontorova 
model with long-range interparticle interactions were studied in [25] . The properties of time pe- 
riodic spatially localized solutions (breathers) on discrete chains in the presence of algebraically 
decaying interactions were described in [32j, [33]. Energy and decay properties of discrete 
breathers in systems with LRI have also been studied in the framework of the Klein-Gordon 
|27j . and discrete nonlinear Schrodinger equations [30] • A main property of the dynamics de- 
scribed by the equation with fractional space derivatives is that the solutions have power-like 
tails. The lattice models with power-like long-range interactions [HI [32], [33], [35], [36], [37J [18] have 
similar properties. As was shown in [TT] [T8] [15] , the analysis of the equations with fractional 
derivatives can provide results for the space asymptotics of their solutions. 

The goal of this paper is to formulate the consistent definition of continuous limit (transform 
operation) for the systems with long-range interactions (LRI). This aim is realized by Proposi- 
tions 1, 4, and Definitions 1, 2. The power-law LRI is considered in [T6] [TT] [181 US] - The exact 
continuous limit results for power-law LRI were formulated in Propositions 2, 3. This operation 
is used to consider a wide class of long-range interactions that can be called alpha-interaction. 
In continuous limit, the equations of motion gives the medium equations with fractional deriva- 
tives. The power-law interaction is a special case of this class of a-interactions. We show how 
the continuous limit for the systems of oscillators with long-range interaction can be described 
by the corresponding fractional equation. 

In Sec. 2, the transform operation that maps the discrete equations into continuous medium 
equation is defined. In Sec. 3, the Fourier series transform of the equations of a system with 
long-range interaction is realized. A wide class of long-range interactions that can give the 
fractional equations in the continuous limit is considered. In Sec. 4, the fractional equations 
are obtained from three-dimensional discrete system. In Sec. 5, the linear power-law long-range 
interactions with positive integer and noninteger powers are considered. The correspondent 
continuous medium equations are discussed. In Sec. 6, the nonlinear long-range interactions 
for the discrete systems are used to derive the Burgers, Korteweg-de Vries, and Boussinesq 
equations and their fractional generalizations in the continuous limit. The conclusion is given 
in Sec. 7. 
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2 Transform operation 



Let us consider a one-dimensional system of interacting oscillators that are described by the 
equations of motion, 

^=gI n ( U )+F(u n ) } (1) 

where s — 1, 2, and u n are displacements from the equilibrium. The terms F(u n ) characterize 
an interaction of the oscillators with the external on-site force. The term I n (u) is defined by 

+oo 

I n {u)= 2j J{n,m) W(u n ,u m ), (2) 

m=—oo 

and it takes into account the interaction of the oscillators in the system. 

For linear long-range interaction we have W(u n , u m ) = u n — u m , and the interaction term 



is 



i n {u) = ^ J{n,m) [u n -u m ]. (3) 



m=—oo 



In this paper, we consider a wide class of interactions (j3J) that create a possibility of presenting 
the continuous medium equations with fractional derivatives. We also discuss the term (j2J) with 
W(u n ,u m ) = f(u n ) — f(u m ) as nonlinear long-range interaction. As the examples, we consider 
f(u) = u 2 and f(u) = u — gu 2 that gives the Burgers, Korteweg-de Vries and Boussinesq 
equations and their fractional generalizations in the continuous limit. 

Let us define the operation, which transforms equations ([1]) for u n (t) into continuous medium 
equation for u(x, t). We assume that u n (t) are Fourier coefficients of some function u(k, t). Then 
we define the field u(k,t) on [-K/2, K/2] as 



(k,t)= ^w^^aKwi, (4) 



where x n = nAx, Ax = 2ir/K is distance between oscillators, and 

i r+K/2 

Un(t) = - dk u(k,t) e te " = J^{u(k,t)}. (5) 

K J -K/2 

These equations are the basis for the Fourier transform, which is obtained by transforming 
from a discrete variable to a continuous one in the limit As — > (K — > oo). The Fourier 
transform can be derived from and (jSJ) in the limit as As — > 0. Replace the discrete 
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u n {t) = {27r/K)u(x n ,t) with continuous u(x,t) while letting x n = nAx = 2im/K — > x. Then 
change the sum to an integral, and equations (TJJ, (jSJ) become 



u(k,t) = dx e- ikx u(x,t) = F{u(x,t)}, (6) 

J — oo 

1 f +co 

u(x, t) = — dk e tkx u(k, t) = ^{uik, t)}. (7) 
2vr > /_ 00 

Here, 

u(k,t) = Cu(k,t), (8) 

and C denotes the passage to the limit Ax — > (K — > oo). Note that u(k,t) is a Fourier 
transform of the field u(x,t), and u(k,t) is a Fourier series transform of u n (t), where we can 
use u n (t) = (2tt / K)u(nAx , t) . The function u(k,t) can be derived from u(k,t) in the limit 
Ax -»■ 0. 

The map of a discrete model into the continuous one can be defined by the transform op- 
eration. 

Definition 1. Transform operation T is a combination T = T~ X C. of the operations: 

1) The Fourier series transform: 

Fa ■ u n (t) -> FA{un(t)} = u(k,t); (9) 

2) The passage to the limit Ax — > 0: 

£ : u(k, t) -> £{m(A;, t)} = u{k, t); (10) 

3) The inverse Fourier transform: 

T~ x : u{k, t) -> ^{u^, t)} = u(x, t). (11) 

The operation T = T~ x £ J 7 a is called a transform operation, since it performs a transform 
of a discrete model of coupled oscillators into the continuous medium model. 

3 From discrete to continuous equation 

Let us consider the interparticle interaction that is described by (jSJ), where J(n, m) satisfies the 
conditions 

oo 

J(n, m) = J{n — m) = J{m — n), l^( n )| 2 < °°- (12) 



n=l 



Note that J(-n) = J(n). 
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Definition 2. The interaction term and ( fH|j m £/ie equation of motion (QJj zs called a- 
interaction if the function 



j a (k)= e ~ lkn J W = 2 E J ( n ) C0S (H ( 13 ) 

satisfies the condition 



n=—oo n=l 



, im [J.(fc) - MQ)] = 



where a > and < IAJ < oo. 



Condition jH]) means that J a (fc) - J Q (0) = 0(|fc| a ), i.e., 

- J«(0) = A,|fc| a + #„(£;), (15) 

for — > 0, where 

lim R a (k)/\k\ a = 0. (16) 

Examples of functions J(n) for a-interactions can be summarized in the table of the Ap- 
pendix. 

Proposition 1. The transform operation T maps the discrete equations of motion 

-^i = g J ^ m ) M*) - + ^M*)) ( 17 ) 

m=— oo 

noninteger a-interaction into the fractional continuous medium equations: 

^u(x, t) - G a A a ^-^u(x, t)-F (u(x, t)) = 0, (18) 

where d a /d\x\ a is the Riesz fractional derivative, and G a = g\Ax\ a is a finite parameter. 

Proof. To derive the equation for the field u(k, t), we multiply equation (ITT)) by exp(— iknAx), 
and summing over n from — oo to +oo. Then 

+ 00 rj s +00 +00 +00 

E e-* knAx —u n (t)=g E e" lfc " A V(n,m) [w n - u m ) + £ e~ iknAx F(u n ). (19) 

n=— oo n=— oo m=— oo n=— oo 

The left-hand side of ( 1191) gives 

y e -iknA X ^Un{t)_ = (f_ ^ -ifcnAx = &u(k,t) 

^ dt s dt s ^ w dt s y ' 



where u(k,t) is defined by (j3J). The second term of the right-hand side of (I19p is 

+00 

]T e -^F( Un )=f A {FW}. 

n=— 00 

The first term on the right-hand side of fTTTTT) is 



rOO +OO 



J] ^ e- lknAx J{n,m)[u n -u r 



n=—oo m=—oc 



+00 +00 +00 +oc 

= ^ e~ iknAx J{n,m)u n - e- lknAx J{n,m)u m . 

n=— 00 m=— 00 n=— 00 m=— 00 

Using (j3J) and (1241) . the first term in r.h.s. of (1221) gives 

+00 +00 +00 +00 

J2 J2 e- iknAx J(n,m)u n = £ e- ifc " Aa; Mn ]T J(m') = t) J a (0), 

n=— 00 m=— 00 n=— oo m'=— oo 



where we use ( I12p and J(m' + n,n) = Jim'), and 



J a (kAx)= e~ iknAx J(n) = jF A {J( n )}. 



n=— 00 



Note that 



+00 +00 +00 +00 

£ ^ e" ifcnA:!; J(n,m)M m = ^ « m ^ e - ifcnAx J(n, m) 

n=~ oo m=—oo m=— oo n=— oo 



hoo +00 



^ M m e- ifcmAa: J] e- ifcn ' Ax J(n') =u(k,t)j a (kAx), 



m=— oo n'=— 00 

n'/O 

where J(m, n' + m) = J(n') is used. 

As a result, equation (fl~9l) has the form 

= 9&(P) ~ UkAx)]u(k,t)+F A {F(u n )}, 
where Ja{^K)} is an operator notation for the Fourier series transform of F{u n ). 



The Fourier series transform T A of (1171) gives (I26p . We will be interested in the limit 
Ax — > 0. Using (EES]), equation (1261) can be written as 



—u(k,t)-G a f aA {k) u(k,t)-F A {F(u n (t))} = 0, (27) 

where we use finite parameter G a = g\Ax\ a , and 

f aA (k) = —A a \k\ a - R a (kAx)\Ax\- a . (28) 
Note that R a satisfies the condition 

R a (kAx) 
hm — r- — : = 0. 

The expression for T aA (k) can be considered as a Fourier transform of the operator ([3]). Note 
that g —>■ oo for the limit Ax — > 0, if G Q is a finite parameter. 

The passage to the limit Ax — > for the third term of (|27|) gives 

£ : J^a^K) -> CF A F(u n ). (29) 

Then 

^a{^K)} = ^{jCFK)} = HF(Cu n )} = F{F(u(x, *))}, (30) 

where we use £JF A = ^-C. 

As a result, equation (1271) in the limit Ax — ► obtains 

— u(A;, - G Q f a (k) u(k, t) - (u(x, t))} = 0, (31) 

where 

u(k,t) = Cu(k,t), f a (k) = CT aA {k) = -A a \k\ a . 
The inverse Fourier transform of (13T|) gives 

^u(x, t) - G a T a (x) u(x, t)-F (u(x, t)) = 0, (32) 



where T a (x) is an operator 



T a (x)=F- 1 {f a (k)} = A a ^. (33) 



Here, we have used the connection between the Riesz fractional derivative and its Fourier 
transform [2]: \k\ a < — > — d a /d\x\ a . 

As a result, we obtain continuous medium equations (fT8l . □ 

Examples of the interaction terms J(n) that give the operators fl33|) in continuous medium 
equations are summarized in the following table. 
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J(n) 


T a (x) 


(^^"<£&M<* + 3/2,l/2, TO) ) 


-d a /d\x\ a 


(-1)™ 

n 2 


-(1/2) d 2 /dx 2 


1 

7T 

n 2 


—in d/dx 


\n\~^ +1 \ (0 < < 2, ^ 1) 


-2T(-P) cos(tt/5/2) dP/dlxf 


l^r (/3+1) , (/? > 2, /3^3,4,...) 


C(P - 1) d 2 /dx 2 


r(l+a/2+n)r(l+a/2-n) > V 2 ) 


-Tiki) da /d\x\ a 


(-1)" 

(1 — Yl 


z sini vrct j ' 


J{n) = l/n\ 


4ei 9/ dx 



8 



4 Fractional three-dimensional lattice equation 

The generalization for the three-dimensional case can be easy realized. Let us consider the 
three-dimensional lattice that is described by the equations of motion 

d s u 

-Qj7~ = 9 J ( n ' m ) [u a -u m ) + F(u a ), (34) 

m=— co 

where n = (ni, n 2 , n 3 ), and J(n, m) = J(n — m) = J(m — n). We suppose that u n (t) are 
Fourier coefficients of the function w(k, £): 

+00 

u(M) = ^ n n (t)e- ikr " = J- A { Mn (t)}, (35) 

n=— 00 

where k = (k\, k 2 , k 3 ), and 



3 

i=l 



Here, are the translational vectors of the lattice. The continuous medium model can be 
derived in the limit |aj| — > 0. 

To derive the equation for w(k, t), we multiply (|34|) by exp(— zkr n ), and summing over n. 
Then, we obtain 



d s u(k,t) _ ^ 



J a (0)-J a (ka) u{k,t)+F A {F(u n )}, (36) 



where J r A{F(u n )} is an operator notation for the Fourier series transform of F(u n ), and 

+00 

J a (ka)= e" ikr "J(n). (37) 

n=— 00 

For the three-dimensional lattice, we define the o-interaction with o = (0:1,02,0:3), as an 
interaction that satisfies the conditions: 

lim [jQ(k) " Ja( ° )] =^, (, = 1,2,3), (38) 

fc^O I 

where < \A ai \ < 00. Conditions (|38|) mean that 

3 3 

J a (k) - J a (0) = ^2A ai \ki\ ai + ^^(k), (39) 

i=l i=l 
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where 

lim i2 aj (k)/|A^| a * = 0. (40) 



In the continuous limit (|a^| — > 0), the a-interaction in the three-dimensional lattice gives the 
continuous medium equations with the derivatives d ai /dx ai , d a2 /dy a2 , and d 03 /dz a3 : 



dS f> t] = - 9 YA jZ^ t] +F(u(r,t)). (41) 

dt s ^ 1 Sih 

1=1 1 1 

This equation describes multifractional properties of continuous medium. 



5 Linear power-law long-range interaction 

Let us consider the chain with linear long-range interaction that is defined by the equation of 
motion 

d s u + °° 

-- g ^ J{n,m) [u n -u m ]+ F(u n ), (42) 



dt s 

m=—oo 



where J(n, m) = J(\n — m\), and 

J(n) = \n\-^ (43) 

with positive integer number (3. 

Proposition 2. The power-law interaction (fJ5D for the odd number (3 is a -interaction with 
a = 1 for (3 = 1, and a = 2 /or /? = 3, 5, 7.... For even numbers (3, (f^) is not a -interaction. 
For odd number (3, the transform operation T maps the equations of motion with the interaction 
( f^gy into the continuous medium equation Ulfy with derivatives of first order for (3 = 1, 

^u(x, t) - iG x |- u(x, t)-F (u(x, t)) = 0, (44) 
and the second order for other odd (3 ((3 = 2m — 1, m = 2, 3, 4, ...), 

—u(x, t)-G 2 — u(x, t)-F (u(x, t)) = 0, (45) 



where 

G\ = ng Ax, G 2 = 1 ' 1 ' B 2m _ 2 g(Ax) 2 (46) 

4(2m — 2)! 



are the finite parameters. 
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Proof. From Eq. (1261) . we get the equation for u(k,t) in the form 

9SU ^ t] + 9 iUkAx) - J a (0)] u(k,t)-F A {F(u n (t))} = 0, 

where 

+oo 

J a {kAx) = e- iknAx \n\' {1+l3) . 

n=—oo 

The function (j4"8l) can be represented by 

UkAx) = J2^tp ( e ~ tknAx + elknAx ) = 2 E cos ( knAx ^ 

n=l n=l 

Then, we can use (section 5.4.2.7 in Ref. [47]) the relations 

cos(nfc) _ (-l) m - 1 (27r) 2m D f fc 

n=l 

where m = 1, 2, 3, and B 2m (z) are the Bernoulli polynomials |46j, which are 



B n (k) = Y,CnB m k r 



m=0 

Here S m are the Bernoulli numbers. Note B 2m -i = for m — 2, 3, 4. 
For /3 = 1, we have 

J a (kAx) - J Q (0) = i(A;Aa;) 2 - ttA;Ax « -tt/cAx. 
For /3 = 2m- 1 (m = 2,3,...), 



Then, 



Jq(A;) = l2m)T (27r)2mjB2m I 2^ } ' ( ° " k ~ 2n) - 



J a (kAx) - J o (0) « 1 f- 1 \, B 2m ^ 2 (kAx) 2 . 

For /3 = 0, we have f[4Tj. section 5.4.2.9.) the relation 

Ecos(nA;) , r , 
i — '- = - In [2 sin(fc/2)] . 

n=l ^ 
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Then, the limit Ax — » gives 

J a (kAx) « - ln(&Ax) -> oo. (56) 

For even numbers /?, 

J a (kAx) - J a (0) /IfcAx^oo (57) 

since the expression has the logarithmic poles. 

The transition to the limit Ax — > in Eq. fj4"T|) with /3 = 1 gives 

9SU ^ t] - Gi fc «(*, - ^{F («(*, «))} = 0, (58) 

where G\ = ugAx is a finite parameter. The inverse Fourier transform of (158]) leads to the 
continuous medium equation (j44j) with coordinate derivative of first order. For s — 1, this 
equation can be considered as the nonlinear Schrodinger equation. 
The limit Ax — > in Eq. (j4"Tj) with /3 = 2m — 1 (m = 2, 3, ...) gives 

+ g 2 fc 2 f) _ T{F t))} = Qj (5Q) 

where G2 is a finite parameter (146j) . The inverse Fourier transform of ( 1591) leads to the partial 
differential equation (T4"5l) of second order. For s = 2, this equation can be considered as a 
nonlinear diffusion equation, and for s = 1 as a nonlinear wave equation. □ 

Proposition 3. T7ie power-law interaction fftfy with noninteger (3 is a-interaction with 
a = (3 for < (3 < 2, and a = 2 for (3 > 2. For < j3 < 2 ((3^1), the transform operation 
T maps the discrete equations with the interaction fcfify into the continuous medium equation 
with fractional Riesz derivatives of order a: 

—u(x,t)-G a A a —^u(x,t) = F{u(x,t)), 0<a<2, (a^l). (60) 

For a > 2 (a 7^ 3,4, 5,.. .), the continuous medium equation has the coordinate derivatives of 
second order 

— u (x,t) + G a ({a-l)—-u{x,t) = F{u{x,t)), a > 2, (a ^ 3, 4, ...). (61) 
ot s o\x\ z 

Proof. For fractional positive a, the function (1481) can be represented by 

+00 

J a (kAx) = J2— ( e ' lknAX + eiknAX ) = Ll i+»( eikAx ) + Li 1+a (^ ikAx ), (62) 



n 1 

n=l 
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where IAp(z) is a polylogarithm function. Using the series representation of the polylogarithm 
Li p {e z ) = r(l - ^{-zf- 1 + Y J C{(3 ~ n) z n , \z\<2n, (3^ 1,2,3..., (63) 



n=0 



we obtain 



J a (kAx) = A a \Ax\ a \k\ a + 2j2 C(1 ~L% 2n) (Ax) 2n (-k 2 r, a / 0, 1, 2, 3..., (64) 



(2n)! 

where ((z) is the Riemann zeta-function, |A;Ax| < 2ir, and 



A a = 2 r(-a) cos f — ) . (65) 



V 2 

From (164")) . we have 

J a (0) = 2C(l + a). 

Then 

J Q (fcAx) - J a (0) = A a \Ax\ a \k\ a + 2j2 ai + *~ 2n \ Ax) 2n (-k 2 r, (66) 

n=l ^ ^ 

where 0,1, 2, 3..., and |A;Ax| < 2n. 
Substitution of (1661) into (|47|) gives 

^^+ fli 4 a |Axri*r«(A ; ,o+ 

+ 2 ^ E C(a | 2 i)T 2n) (^) 2 "(-fe 2 )"^(^ *) - -^a{F KW)} = 0. (67) 
We will be interested in the limit Ax — > 0. Then, Eq. fl67|) can be written as 

-u(k,t) + G a T aA (k)u(k,t)-F A {F(u n (t))} = 0, 0,1,2,..., (68) 

where we use the finite parameter 

G a = g\Ax\ min{a ' 2} , (69) 



and 



A a |fc|--|Ax| 2 - a C(«-l)^ 2 , 0<a<2, (a^l); 
"" A ^ |Aa;| Q -M a |fc| a -C(a-l)fc 2 , a > 2, (a ^3,4,...). 1 J 



The expression for 7a, a(&) can be considered as a Fourier transform of the interaction operator 
([2]). From (169]) . we see that g — > oo for the limit As — > 0, and finite value of G a . 



13 



The transition to the limit Ax — > in Eq. f)68p gives 

—u(k, t) + G a f a (k)u(k, t) - F{F (u(x, t))} = (a + 0, 1, 2, ...), (71) 



where 



T(k) 1 Aa \ k \ a > 0<a<2, a ^ 1; 

" ' S ' 1 <(a- 1) A; 2 , 2 < a, a ^ 3,4,.... 1 ' 



The inverse Fourier transform to ( 17TT) is 
9 



u(x,t) + G a T a (x) u(x,t) -F(u(x,t)) = a ^0,1,2,..., (73) 



where 

' -A a d a /d\x\ a , (0 < a < 2, a / 1) 
C(a-l) <9 2 /9|x| 2 , (a>2, a ^3, 4,.. 

As the result, we obtain the continuous medium equations (1601) and (|6TT) . □ 



r a (x) = ^- 1 {t(fc)} 



For s = 1 and -F(u) = 0, Eq. fl60l) is the fractional kinetic equation that describes the 
fractional superdiffusion [121 US], [39] . If F{u) is a sum of linear and cubic terms, then Eq. fl60l) 
has the form of the fractional Ginzburg-Landau equation [40l I4T1 1421 l43l 144] . A remarkable 
property of the dynamics described by the equation with fractional space derivatives is that 
the solutions have power-like tails. 



6 Nonlinear long-range interaction 

Let us consider the discrete equations with nonlinear long-range interaction: 

In{u)= ^2 Ja(n,m)[f(u n )-f(u m )], (74) 

m=— oo 

where f(u) is a nonlinear function of u n (t), and J Q (n, m) = J a (n — m) defines the a-interaction. 
As an example of J Q (n), we can use 

(-1)" 

Ja{n) = T(l + a/2 + n)T(l + a/2-nY (?5) 

For a = 1,2,3,4, the interactions with f{u) = u 2 and f{u) = u — gv? give the Burgers, 
Korteweg-de Vries and Boussinesq equations in the continuous limit. For fractional a in Eq. 
r5|) . we can obtain the fractional generalization of these equations. 
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Proposition 4. The transform operation maps the equations of motion 



d s u n {t) 



+00 



9 ^ J a( n - m )[f( U n) - f(Um)]+ F(u n ), (76) 



dt s 

m=—oo 

where F is an external on-site force, and J a {n) defines the a -interaction, into the continuous 
medium equations 

where G a = g\ Ax\ a is a finite parameter. 

Proof. The Fourier series transform of the interaction term (174"|) can be represented as 



-oo +00 



e- iknAx I n (u)= Yl E e- iknA *J(n,m)[f(u n )-f(u m )} 



n=—oo m=— 00 



+00 +00 +00 +00 



E E e- ik " A *J(n,m)f(u n )- ^ £ e" toA V(n, m)f(u m ). (78) 



n=— oo m=— oo n=— oo m=—oo 

For the first term on the r.h.s. of (1781) : 

+00 +00 +00 +00 



E E e-* knA *J(n,m)f(u n )= ^ e~ iknAx f{u n ) ^ J(m>) = F A {f(u n )} J a (0), 

(79) 



n=— 00 m=—oo n=—oo m =—oo 



where we use J{m' + n,n) = J{m'). 

For the second term on the r.h.s. of (1781) : 

+00 +00 +00 +00 

E E e-* knA *J(n,m)f(u m )= £ /(« m ) £ e -^J(n,m) = 

n=—oo m=—oo m=—oo n=—oo 

+00 +00 

= E f( u ™) e ~ tkmAx E e" ifcn ' A V(n') = ^a{/K)} L{kAx), (80) 



n =— 00 



where we use J(m, n' + m) = J(n'). 
As a result, we obtain 



= g[J a ( Q ) _ J Q (^)]^ A {/K)} + f A {FW}, 
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where u(k,t) = J r A{u n (t)}, and J a (kAx) = J- A {J{n)}. 
For the limit Aa; — > 0, Eq. flHTj) can be written as 



—u(k,t)-G a f aA (k) u(k,t)-F A {F(u n (t))} = 0, (82) 

where we use finite parameter G a = g\Ax\ a , and 

T aA (k) = -A a \k\ a - R a (kAx)\Ax\~ a . (83) 
Here, the function R a satisfies the condition 

R a (kAx) 
hm — — — : = 0. 

Az^O \Ax\ a 

In the limit Ax — » 0, we get 

— t) - G Q f a (k) T{f (u(x, t))} - F{F {u(x, t))} = 0, (84) 

where 

u(k,t) = Cu(k,t), TJJz) = Cf aA (k) = -A a \k\ a . 
The inverse Fourier transform of fl84"j) gives the continuous medium equation (I77|) . □ 

Let us consider examples of quadratic-nonlinear long-range interactions. 

1) The continuous limit of the lattice equations 

du (t) +0 ° +0 ° 

Q f = 9i J 1 (n,rn)[ul-u 2 m ]+ g 2 J 2 (n, m)[u n - u m ], (85) 

m=— oo m=—oo 

where Ji(n) (i = 1,2) define the ^-interactions with a\ — 1 and ct2 = 2, gives the Burgers 
equation [38] that is a nonlinear partial differential equation of second order: 

d d d 2 

— u(x, t) + Giu(x, t)—u(x, t) - G 2 -^u(x, t) = 0. (86) 

It is used in fluid dynamics as a simplified model for turbulence, boundary layer behavior, shock 
wave formation, and mass transport. If we consider J 2 (n,m) with fractional a 2 = a, then we 
get the fractional Burgers equation that is suggested in [19]. In general, the fractional Burgers 
equation is 

d d ai d a2 

— u(x, t) + G ai u(x, t)— — -u(x, t) - G a2 -^u(x, t) = 0. (87) 
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2) The continuous limit of the system of equations 
du n (t) 



91 Yl J l( n ' m )t M n - "ml + 93 Yl J 3( n i m )l u 



dt 

m=—oo m=—oc 

where Ji{n) {i = 1,3) define the aj-interactions with a>i = 1 and = 3, gives Korteweg-de 
Vries (KdV) equation 

^u(x, t) - G lU (x, t)-^u(x, t) + G 3 -^u(x, t) = 0. (89) 

First formulated as a part of an analysis of shallow-water waves in canals, it has subsequently 
been found to be involved in a wide range of physics phenomena, especially those exhibit- 
ing shock waves, travelling waves and solitons. Certain theoretical physics phenomena in the 
quantum mechanics domain are explained by means of a KdV model. It is used in fluid dynam- 
ics, aerodynamics, and continuum mechanics as a model for shock wave formation, solitons, 
turbulence, boundary layer behavior and mass transport. 

If we use noninteger aj-interactions for Ji(n), then we get the fractional generalization of 
KdV equation [501 EI]: 

d d ai d a3 

— u(x, t) - G ai u(x, t)-^-^u(x, t) + G a3 -g-^u(x, t) = 0. (90) 
3) The continuous limit of the equations 

d 2 u n (t) +oc + °° 

o" 2 = 92 Y M n ' m )lf( U n) - f(Um)] +94 ^ M n i m )l U }, (91) 

m=— oo m=—oo 

where 

f( U ) = U - g U 2 , 

and Jiin) define the cv interactions with a 2 = 2 and a 4 = 4, gives the Boussinesq equation 
that is a nonlinear partial differential equation of forth order 

d 2 d 2 d 2 <9 4 

0pu{x, t) - G 2 -^u(x, t) + gG 2 -^u 2 (x, t) + G 4 -^u{x, t) = 0. (92) 

This equation was formulated as a part of an analysis of long waves in shallow water. It was 
subsequently applied to problems in the percolation of water in porous subsurface strata. It also 
crops up in the analysis of many other physical processes. The fractional Boussinesq equation 
is 

Q2 QOL2 QCX2 QOti 

Qpu{x,t) - G a2 —u(x,t)+gG a2 —u 2 (x,t) + G a4 —u(x,t) = 0. (93) 



17 



7 Conclusion 



Discrete system of long-range interacting oscillators serve as a model for numerous applications 
in physics, chemistry, biology, etc. Long-range interactions are important type of interac- 
tions for complex media. An interesting situation arises when we consider a wide class of 
a-interactions, where a is noninteger. A remarkable feature of these interactions is the exis- 
tence of a transform operation that replaces the set of coupled individual oscillator equations 
by the continuous medium equation with the space derivative of noninteger order a. Such a 
transform operation is an approximation that appears in the continuous limit. This limit allows 
us to consider different models in unified way by applying tools of fractional calculus. 

We can assume that an asymmetric interaction term (J(n — m) 7^ J(\n — m\)) leads to other 
forms of the fractional derivative [2]. 

Note that a fractional derivative can be results from a fractional difference as interaction 
term, just as an nth order difference leads to an nth derivative [2]. It follows from the repre- 
sentation of the Riesz fractional derivative by Grunwald-Letnikov fractional derivative: 

d a , v 1 1 A (-lVT(a + l) . . , . , .. 

- , uix.t) ~ : r^~r- > — ; ; — ; r \u(x — nh, t + u [x + nh, t , 

d\x\ a V ; 2cos(7ra/2) /i" ^ T(n + 1) T(a-n + l) 1 V ; V n 

(94) 

where h = Ax is the discretization parameter. 

A similar approach to fractional dynamics in the context of the diffusion equation was 
developed in the papers [52j [53]. In those papers, a continuum limit of (non- interacting) 
random particle motions leads to a fractional equation. 
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